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Abstract 



The theory of sound propagation in macroscopically isotropic and 
homogeneous porous media saturated with supcrfluid ''He has been 
developed neglecting all damping processes. The case when the normal 
Cu I fluid component is locked inside a porous medium by viscous forces is 

investigated in detail. It is shown that in this case one shear wave and 
two longitudinal, fast and slow, waves exist. Fast wave as well as slow 
wave is accompanied with temperature oscillations. The velocities of 



Q I these waves are obtained. 



PACS numbers: 81.05.Rm; 47.35.+i; 82.33.Ln; 47.37.-hq; 
KEY words: porous materials; hydrodynamical waves; 



1 Introduction 

The investigation of acoustic propagation in a fluid-saturated porous me- 
dia provides opportunities to study disordered pore structures, as well as 
the properties of the fluids confined in the pores. In [1] Blot developed a 
simple theory of propagation of elastic waves in statistically homogeneous 
and isotropic porous media saturated with a (one-component) compressional 
fluid assuming that pore sizes are smaller compared to the wavelength i.e. 
the long wavelength limit is considered. Blot's theory predicts one shear 
wave and two longitudinal, fast and slow compressional waves. The fast 
compressional wave corresponds to solid and fluid moving in phase, and the 
slow one corresponds to solid and fluid moving in opposite phase. The Blot's 
theory was experimentally confirmed by the observation of the fast and slow 
compressional wave in glass bead samples filled with water [2]. 

Because of the unique properties of liquid helium [3], it is of funda- 
mental importance to investigate the sound propagation in superfluid-filled 
porous media. Below the A-transition temperature T\ = 2.17i^, liquid *^He 
becomes a superfluid and behaves as a "mixture" of a normal(viscous) com- 
ponent with mass density pn and a superfluid (zero-viscosity) component 
with mass density ps', it is clear that Pn + Ps = Ph where ph is the mass 
density of superfluid *^He(so-called He II). It must be pointed out that the 
superfluid component and the normal component can not be physically sep- 
arated [3]. The normal (superfluid) component fraction increases (decreases) 
monotonously from p„ = 0(/9s = P/i) at T = to /9„ = ph{ps = 0) at and 
above T\. It is important that p„ and ps can move independently, and 
therefore in bulk He II the two mode of sound propagation are possible. 
In first sound the two fluid components oscillate in phase; and in second 
sound ps and p„ oscillate in opposite phase. One of the most interesting 
property of He II is fourth sound mode. In fourth sound only the superfluid 
component moves, when the normal fluid component is locked inside a rigid 
porous medium by viscous forces, so that its velocity is identically zero. 

Since a superfluid-filled porous medium is really a three-component sys- 
tem. Blot's theory is not applicable directly to this system. But, at low 
temperatures (T < l.\K), when the normal fluid component fraction may 
be neglected, Johnson has shown that fourth sound mode is exactly a Blot's 
slow wave [4]. It should be noted that at low temperatures, when the mean 
free path of phonons and rotons becomes larger than the pore size the hy- 
drodynamical description breaks down, so the comparison to Biot's model 
should be made at higher temperatures [5]. 



We have to point out that a number of authors extended Biot's theory in 
the case of two or more fluids (for review see [6] , [7] ) . Since the fluids usually 
considered are not superfluids, the authors neglected some effects, which are 
actually important in the case of superfluids. Indeed, it is well known that 
in the case of ordinary fluids the fountain(and/or mechano-caloric) effect is 
negligible, however, in the case of a superfluid these effects are crucial. In 
particular, a temperature gradient produces relative motion between the two 
fluids in superfluids [3]. Therefore we have to distinguish between pressures 
with a "mechanical origin" and those with a "thermal origin" , the later can 
only be found in superfluids [3] whereas the former can appear in any fluid. 
Moreover, it is important that the superfluid component and the normal 
component can not be physically separated. Actually He II is the fluid with 
additional degrees of freedom [3], but not a real mixture of two fluids. 

In this paper the generalization of Biot's theory is given in the case 
of superfluid-saturated porous media, when all damping processes are ne- 
glected. The organization of this paper is as follows. In Section II we apply 
Biot formalism to bulk He II and derive some equations, which are used 
in Section IV. Section HI and Section IV give a derivation of the equa- 
tions which describe the propagation of longitudinal and shear waves in a 
superfluid-saturated, homogeneous and isotropic porous media. This is the 
main result of this paper. In Section V the formalism is applied to the most 
important particular case when the normal fluid component is locked inside 
a porous media by viscous forces. Finally, in Section VI we discuss some 
limits of our equations and summarize results. 

2 The case of bulk He II 

In this section we apply formalism developed by Biot [1],[8] to the wave pro- 
cesses in bulk He II. The wave processes in bulk He II is well understood [3] , 
but Biot's formalism has the advantage that it can be generalized easily for 
three-component systems. Besides, in this section we derive some equations 
which are used later. 

The basic idea is to describe the wave processes in He II in terms of 
"elastic deformations" of the superfluid and the normal components of He II. 
In this formalism the pressure with a "thermal origin" have to be considered 
as an "elastic" restore force as well as the pressure with a "mechanical 
origin". Here we present the key points of calculations, for more discussion 
of details of the formalism the reader is referred to [1],[8]. 



Let us define the displacement of the superfluid and the normal com- 
ponents as Us and Un respectively. Then the displacement u of He II is 

[3]: 

U = Us H Un- (1) 

Ph Pn 

Following [1] we choose the Lagrangian coordinates of the system as the 
six displacement components Usx,Usy,Usz,Unx,Uny,Unz of He II. Note that 
in contrast with ref.[l] no averaging procedure is required here. Then the 
kinetic energy Th of the unit volume of He II is presented as: 



2T/J = ^ Pn 






x,y,z \ 

here we take into account that the superfluid component and the normal 
component move independently. So the mass coupling parameter between 
the normal component and the superfluid component equals zero [1]. 

Let us assume that the system is in equilibrium being at rest. Therefore, 
any displacement is a deviation from a state of minimum potential energy. 
Then the potential energy W^ per unit volume of He II is [1],[9]: 

2Wh = a^en + (T%; (3) 

here we take into account that the shear modulus of He II is zero. We 
formally separate the stress of the normal component cj„ and the stress of 
the superfluid component a^ (see Section III as well); they describe the 
pressure with a "thermal origin" as well as the pressure with a "mechanical 
origin", e^ = div Ug and e„ = div Un are the strain in the normal and the 
superfluid components of He II respectively. From (1) the strain in the fluid 
e is given by: 

Ps , Pn /.N 

e = —es H e„. (4) 

Ph Ph 

When there is no (nonuniform) relative motion between the superfluid 
and the normal components then e^ = e„ and no temperature gradients exist 
(see the text after equation (16)). In this case the dynamics of the system 
is described by the theory of elasticity. 

But when e^ 7^ e„ temperature gradients exist and the fountain and/or 
mechano-caloric effect takes place (see the text after equation (16)). We 
show in this section that in this case the formalism of theory of elasticity is 
applicable as well. 



In the first approximation the stress components are linear functions of 
the strain components, so one gets [1]: 

Ug — n,gts -r J^gntnj 

'^n — ^n^n + Rgn^s] (5) 

Rg and R[[ describe "elastic properties" of the superfluid and the normal 
components, and R^^ describes the "elastic interaction" between them [1]. 

To make results more useful let us express Rg, i?„ and i?g„ in terms of 
directly measurable coefficients [8]. 

Let us now consider a volume of He II in a vessel. First we consider the 
case, when none of the fluid components is allowed to escape from vessel. We 
shall refer to this case as "unjacketed compressibility test" (see Section III). 
If a pressure p is applied to He II (see Fig. 1), no relative motion between 
the superfluid component and the normal component takes place, i.e.: 

eg = en = e. (6) 

So, in this case temperature difference does not exist. Fluid pressure has 
only "mechanical origin" and we have [8]: 

h Ps 

cTs = — p; 



Ph 

Pn 

1 

Ph 



ol = -^p; (7) 



ii.g€g -]- Kg^en', 



so, from (5) and (7) we get: 

Ps 

P 

Ph 

P = Rn^n + Rsn^s] (8) 

Ph 

Besides, in this case we have: 

1 _ _e 

Kf p' 

where Kf is the adiabatic bulk modulus of He II. 

Using (9) from (8) we get (recall that e = e„, in this case) 



(9) 



yh _ Ps_T^^ jyh 

Ph 



R, — — Kf - Rgj^. 



Ri = ^Kf-Rl. (10) 

Ph 



Now we consider another case, when the superfluid component is allowed 
to escape through a superleak, but not the normal fluid component (see 
Fig. 2). We shall refer to this case as "second jacketed compressibility 
test" (see Section III). If a pressure p is applied to He II, the superfluid 
component passes through the superleak into the reservoir of He II. An 
applied external pressure squeeze out some of the superfluid component, this 
means that temperature increases in a vessel. So, a temperature difference 
between the vessel and the reservoir of He II exists. As it is well known, 
temperature difference between the two volumes of He II is accompanied 
by a pressure difference between these volumes (fountain or/and mechano- 
caloric effect) [3]. Since the superfluid component is allowed to pass through 
the superleak, its partial pressure is zero. But the normal fluid component is 
identified with the thermal excitations of He II (gas of phonons and rotons) 
and its partial pressure is not zero. Thus, if we neglect the thermal expansion 
coefficient of He II (which is anomalously small) , we can conclude that fluid 
pressure in this case has only "thermal origin". It is clear that the normal 
fluid component pressure equals to the external pressure (let us define it in 
this case by p'), and we have: 

^s = 0; 

^:^ = -P'; (11) 



or using (5) and (11): 



— Rq^s ~l" Rtsn^ 



-p' 



R-n^n + R«r,^s\ (12) 



Let US define the bulk modulus of the normal fluid component Kn as: 

^ = -^; (13) 

Kn V 

Since the superfluid component and the normal component can not be phys- 
ically separated, the deflnition of the bulk modulus of the normal fluid com- 
ponent is formal. But let us point out that K^ does not enter into final 
equations (18), (19); so this problem is irrelevant in this case. 
From (13) and (12) we have: 

{R\S^R\Kn{\-^=^. (14) 



Besides, in this case the energy conservation law gives: 

-^{en-es)SphT = phC^T- (15) 

Ph 

where {ps/ Ph){^n — ^s) represents the fluid volume entering through the 
superleak [8]. C and S are the specific heat and the entropy per unit mass 
of He II, respectively. AT is the difference between the temperature of the 
fluid inside a vessel and the temperature of the reservoir of He II caused by 
external pressure [3]: 

^^ = ^- <'<=) 

One can see from (15) that AT = when e^ — e„ = 0, so pressure has only 
a "mechanical origin". But when e^ ^ e„ from (15) one gets that AT ^ 0, 
therefore temperature gradients exist. 

From (15) using the first equation in (12) we can eliminate e^ and e„, 
and we get the following expression for i^„: 

here we have used (10), (14), (16) and definition (13) as well. 
From (10) and (14) using (17) we directly get: 

PhS^Tpl 



R's - 


Ph 


T>h 


= ""^ Pi 


^sn = 


Ph 



+ 



+ 



c pV 

PhS^T pl 

c pV 

PhS^T pl 

c pI 



(18) 



Equation (18) gives Rg, i?„ and Rg^ in terms of directly measurable 
coefficients. These relations are used in Section IV. 

Now, to demonstrate the correctness of the expressions obtained above 
let us calculate the velocity of sound propagation in He II. 

Using (2), (3) and (5) from Lagrange's equations one derives the equa- 
tions for the wave propagation in He II [1]: 

q2 
Ps-Q^^s = V (i?^es + i?^„e„); 

Pn-g^en = W\R':^es + Rten); (19) 



Solutions of (19) with k wave vector have the form: 

Cs = Csex-p[ik{f+Vt)]; 

en = Cnexp[ik{f+Vt)]; (20) 

Cs and Cn are real constants, and V is the velocity of the wave. 

Substituting (20) into (19) we find that there are two modes of longitudi- 
nal waves the fast and the slow waves with velocities Vi and V2 respectively. 
Using (18) we get the expressions for Vi and V2: 



V^' - 


5 




Ph 


Vi - 


PsS^T 
n 



(21) 

The fast wave corresponds to the superfluid component and the normal 
component moving in phase; from (19) one can see that in the case of fast 
wave Cs = Cn [1]- This means that e^ = e„, so no fountain(and/or 
mechano-caloric) effect takes place and only the pressure with a "mechanical 
origin" appear. This is a ordinary wave in an ordinary fluid (see also Section 

V). 

The slow wave corresponds to the superfluid and the normal components 
moving in opposite phase. In the case of slow wave (19) gives that psCg + 
PnCn = [1]. Thus e = and pressure with a "mechanical origin" does 
not appear; but e^ 7^ e„ , thus temperature oscillations exist. This is a 
temperature wave which appear only in superfluids (see also Section V). 

These are well known results of Landau's theory of sound propagation in 
He II [3]. Therefore, we have shown that the formalism developed by Biot 
[1] may be applied to the wave processes in He II. Besides, we conclude that 
bulk He II may be regarded as a fluid-saturated porous media with porosity 
ps/ Ph- The normal fluid component plays the role of the solid part with Kn 
bulk modulus and zero shear modulus, and the superfluid component plays 
the role of the fluid. Equation (21) shows that second sound is exactly slow 
wave mode. 

3 Derivation of general equations 

In this section we will derive the general equations which govern the propaga- 
tion of longitudinal and shear waves in a superfluid-saturated, homogeneous 
and isotropic porous media. 



Let us consider a unit cube of superfluid-filled, macroscopically homo- 
geneous and isotropic porous solid as an element. The element is assumed 
to be large as compared to pore sizes, but small compared to a wavelength. 
In this case, the element is described by the average displacement of the 
normal fluid component ■«„, the superfluid component -Ug, and the solid u 
parts. 

Following [1], let us introduce the Lagrangian viewpoint and the concept 
of generalized coordinates. The Lagrangian coordinates are chosen as the 
nine average displacement components of both the solid and the fluid i.e. 

The kinetic energy T of the system per unit volume is [1]: 



2r 



E 

i=x,y,z 
2P?2 



Pn 



duj 
dt 



dui dus 



dt dt 



+ P22 

+ 2/9?2 



dt 



+ P22 



dUr, 



dt 



+ 



dUj dUni 

dt dt 



(22) 



here we take into account that the directions x, y, z are equivalent and un- 
coupled dynamically. We can express pn, p|2' P22' P12 ^"^^ P12 ^^ terms 
of directly measurable coefficients. Calculations quite similar to [1] for pn, 
^22' P221 Pi2^ P12 density coefficients give (see also [4]): 



Pn = 


= {l-(3)psol- 


-PI2- 


-PI2. 


P22 = 


= (^Ps-Pu, 






P^22 - 


= Ppn-Pu; 






P12 -- 


= {l-a)pps; 






n 

P12 = 


= (1 - a)/3p„; 







(23) 



From (23) one can calculate pn, P22, P22) P12 ^^'^ P12 using the directly 
measurable coefficients. /3 is the porosity(fluid volume fraction); p^oi is the 
density of the solid; a is a purely geometrical quantity independent of solid 
and fluid densities; a = r? where n is the index of refraction [4]. 

The Pi2{Pi2) coefficient represents a mass coupling parameter between 
the solid and the superfluid (normal) components. As ps and p„ can move 
independently, the mass coupling parameter between the superfluid and the 
normal components is identically zero. 

Let us assume that system is in equilibrium being at rest, therefore any 
displacement is a deviation from a state of minimum potential energy. Then 




the potential energy W per unit volume of the system is given by [1]: 

2W = axCx + ayCy + U^e^ + TxJx + Ty-fy + Tzjz + CTses + (Tnen, (24) 

here, the generalized procedure of the classical theory of elasticity has been 
used [9]. W is expressed in terms of the strain and stress components. 
The stress tensor in the porous material is: 

(^x + CTs + O-n Tz 

Tz ay + as+ (Jn 

ax,cry, az,Tx,Ty,Tz components of the stress tensor are the forces applied to 
the portion of the cube faces occupied by the solid. Note, that the fraction 
of solid area per unit cross section is (1 — /3) [1]. In stress tensor we formally 
separate the forces ag and (T„ acting respectively on the normal component 
and the superfluid component "parts" of each face of the cube. The "frac- 
tion" of the superfluid(normal) component area per unit cross section is 
PPs/Ph iPpn/Ph)- Besides, a = as + (T„; a is the force acting on the fluid 
part of each face of the cube. The fraction of the fluid area per unit cross 
section is /3. 

The strain tensor in the solid is: 

Sx 2lz 2^y 
olz Hy 2^z 

with ex = dux/dx etc. and ^z = dux/dy + duy/dx etc. 

In the first approximation the stress components are linear functions 
of the strain components. The stress-strain relations for superfluid-filled, 
statistically isotropic porous media can be expressed as [1]: 

ax = 2Nex + Ae + QsCg + Qn^^n] 
ay = 2Ney + Ae + QgCg + Qnen, 
az = 2Nez + Ae + QgCg + Qnen, 

Tx = N-fx] Ty = N-fy] Tz = Njz] 



where 



^x ~r Gy -p 62 

10 



Note that, because of (25) describes elastic properties of coupled system 
of the porous solid and He II, the coefficients Rg, Rn, Rsn differ from R^, 
i?^, Rg^. But the relations derived in Section II are still satisfied. Indeed, 
it is clear that one can consider porous solid as a vessel. So, we can iden- 
tify "unjacketed compressibility test" and "second compressibility test" in 
Section II with the unjacketed compressibility test and the second jacketed 
compressibility test considered in the next section respectively. In Section 
IV we use some of the relations derived in Section II. 

It should be pointed out that thermodynamic parameters of He II in the 
porous media and the ones of bulk He II may differ from each other (/?„ for 
example) [3]. 

The thermal expansion coefficient of He II is too small and we can neglect 
it. Besides, we neglect the heat exchange between the solid and He II. This 
is a reasonable assumption if the wave frequency is not too small. Under 
this assumptions it is clear that between the solid and He II forces caused 
by temperature gradients do not exist. So, forces between the solid and He 
II have "mechanical origin". The two last terms in the first equation from 
(25), for instance, describe the interaction of the superfluid component and 
the normal component with the solid [1]. Because of He II interacts with 
the porous solid "as a whole" we have: 

and using (4) one gets: 

Qs = —Q; Qn = —Q; (26) 

Ph Ph 

Now we are ready to write down the general wave equations. The cal- 
culations following the procedure of ref. [1] gives the equations for the wave 
propagation in a superfluid-saturated porous media. Using (22), (25) and 
(26) from Lagrange's equations we derive: 

NV'^u + {A + N)Ve + ^QVe, + —QVen = 

Ph Ph 

= -Q^{PW^ + Pl2'^s + PuUn); 

—QVe + RsVes + Rsn^en = wr^{Pi2^ + P22^s)] 

Ph Ot'^ 

^QVe + RnVen + RsnVes = |^(p?2^+ P22^n); (27) 

Ph ot^ 

11 



here V is the nabla operator. 

Because of the statistical isotropy of the material, the rotational waves 
are separated from the longitudinal waves and obey independent equations 
of propagation [1]. 

Applying the divergence operator to the set of equations (27) we obtain: 



2 ^ D. , /'^ ^. , /'^ 



92 



V'[Pe + ^Qes + —Qen] = ^{pne + p'^es + Pu^n); 
72 fP- 



Ph Ph ) oi^ 

92 



V — Qe + iJses + i?s„e„ = -Tr^{p\2e- ^ P22^s)\ 



V2 f ^Qe + i?,„e, + i?„e„) = ^K2e + P22en); (28) 



with 



Ph ) dt^ 

Qe + Rsnes + Rnen] = TTT^i 

Ph J dV 

P = A + 2N 



Equations (28) govern the propagation of longitudinal waves. 
Applying the rot operator to equations (27) we obtain: 

Q2 

NV^UJ = —2 {puUJ + pl2'^s + Pu'^n); 

= Q^iPl2^ + P22^s); 

= ^{p^^u + p-^Con); (29) 

with the definitions 

uj = rot u; uJs = rot Ug] COn = rot Un, 

Equations (29) govern the propagation of pure rotational(shear) waves. 

When ps = or pn = equations (28), (29) reduces to the set of equa- 
tions obtained by Biot [1]. Besides when (3 = 1 and Kb = N = (the case 
of bulk He II) (28), (29) reduce to (19). 

Equations (28), (29) describe the propagation of longitudinal and shear 
waves in a superfluid-saturated, homogeneous and isotropic porous media 
in the long wavelength limit when all damping processes are neglected. 

Note that using results the of [1] and [5] one can easily incorporate 
the effect of the wave attenuation due to finite viscousity of the normal 
component in these equations. 



12 



4 Calculation of the coefficients 

To make the theory useful we must express coefficients A'', A, Q, Rg, Rn, Rsn 
in terms of directly measurable ones [8] . 

The coefficient N represents the shear modulus of the bulk material, 
which is directly measurable. A, Q, Rg, Rn, R^n can be expressed in terms of 
directly measurable coefficients by considering three compressibility tests. 
In the absence of electrochemical interfacial effect between fluid and solid 
A,Q,Rsn,Rs,Rn and A^ are independent of what fluid is in the pores, in- 
cluding vacuum [4]; we assumed this through this paper. 

In the "unjacketed" compressibility test, a sample of the material is 
immersed in a fluid to which a pressure p is applied (see Fig. 3). Then the 
adiabatic bulk modulus of solid Ks is determined by: 

1 e 

In this case (6) takes place, thus fluid pressure has only "mechanical 
origin" and in (25) we have [8]: 

f^x = (yy = (Tz = -(1 - 0)p] 

oPs 

o-g = -p—p; 

Ph 

CTn = -P^P. (31) 

Ph 

And from (25), (26) and (31) we get: 

-{l-P)p = lNe + Ae + ^Qes + —Qen; 
3 Ph Ph 

-P—p = —Qe + RgCs + Rsnen, 
Ph Ph 

-P—p = —Qe + Rsn^s + Rnf-n, (32) 

Ph Ph 



Using (30) from (32) we get: 



{A + {2I?,)N)Z,-Q^Z2 
""' = {l-P)Z,-PQZ, ' ('') 



where 



sn) ) 



Z\ — RgRn — {R. 

p p p p 

Z2 = —nRn — 2 n-Rsn H ^-Rs] (34) 

Ph Ph Ph 
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Besides, (30) and the two last equations from (32) give: 



'snt 



V KsJ Ph 

Rn = (p-^)—Kf-Rsn; (35) 

V J^S/ Ph 

here we have used (6) and (9) as well. 

In the first "jacketed" compressibility test, a sample of the material is 
enclosed in a thin impermeable jacket and then subjected to an external 
fluid pressure p. The fluid inside the jacket can escape through a tube in 
the reservoir of He II (see Fig. 4), so the fluid pressure is exactly zero. The 
adiabatic bulk modulus of the skeletal frame Kb is: 

1 _ e 
Kb p' 

and in (25) we have [8]: 



(36) 









O'x = 


ay = a, = - 


-p; 








(Js = 


CT„ = 0. 




And from 


(25), (26) 


and 


(37) we 1 


lave: 








p -- 


= ^Ne + Ae + ^Qe, 
3 Ph 


5+ Qe 
Ph 






= 


Ph 


+ Rs^-s + Rsr 


i^nj 






= 


- f^Qe 
Ph 


+ Rsn^s + Ri 


i^n) 


Using 


(36) from 


(38) 


we get: 










Ki 


(^ + 


■ (2/3)iV)Zi - 


Q'Z2 



(37) 



(38) 



(39) 

where Zi and Z2 are given by (34). 

In the second "jacketed" compressibility test only the superfluid compo- 
nent is allowed to escape through a superleak in the reservoir of He II (see 
Fig. 5). An applied external pressure squeeze out some of the superfluid 
component; so the temperature of He II in the pores increases and mechano- 
caloric effect [3] takes place. In this case the bulk modulus of the skeletal 
frame with the normal fluid component Kh is defined by: 

7^ = -'-^ (40) 

Kh P 

14 



Since only the superfluid component is allowed to escape through the super- 
leak, the effective frame here includes both the solid and the normal fluid 
component. Because of the superfluid component and the normal compo- 
nent can not be physically separated, the definition of the bulk modulus of 
the effective frame is formal. But let us point out that K}j does not enter 
into final equations (44); so this problem is irrelevant in this case. 

Let us note again that we neglect the heat exchange between the porous 
solid and He II. In general this is not the case. But recall that we consider 
the wave processes in a porous medium. And, if the wave frequency is not 
too small, the heat exchange between the porous solid and He II during 
times of order of the wave period is negligible (adiabatic process). Besides, 
we neglect the thermal expansion coefficient of He II (see Section II). So, 
the response of the sample on the external pressure p simply is a sum of the 
purely "mechanical response" of the solid and the purely "thermal response" 
of the normal fiuid component. Thus, we can represent the external pressure 
as follows: p = {p — fip') + (3p' . The first term in this sum corresponds to the 
purely "mechanical response" of the solid, and the second term in this sum 
corresponds to the purely "thermal response" of the normal fiuid component. 
In this case in (25) we have: 

cfz = (Jy = az = -{l - (3')p; 
as = 0; 

-P'p; (41) 



O". 



n 



here /?' = P{p' /p) is the measure of the "thermal response" of the effective 
frame. If we neglect the mechano-caloric effect /3' = and this case reduces 
to the first compressibility test. 
From (25), (26) and (41) we get: 

-(1-/3')P = lNe + Ae + ^Qes + —Qen; 
3 Ph Ph 

Ps 
= Qe + Rs€s + Rsnen, 

Ph 

-f3'p = ^Qe + Rsnes + Rnen, (42) 

ph 

Using (40) and two last equations from (42) we can write: 

n - P'n P ^R ^^" P' p P' ■ 

Ph Kh R'^ Kn Kn 
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here we have used (13) and (12) as well. 

Finally (33), (39) with (35), (43) and (17), (18) directly give: 

^ ^ (1 - /?) (1 - /? - {Kb/Ks)) + (5{Ks/Kf)KB 2^^^ 



Q 



l-P- {Kb/Ks) + ms/Kf) 3 

{l-(i-{KB/Ks))(3Ks 



I -13- (Kb/Ks) + PiKs/Kf) ' 



I- P- {Kb/Ks) + ms/Kf) pI^^ C pf 

„ ^ flKs_ pI PhS^T pl 

I- P- {Kb/Ks) + ms/Kf) pI^'' C pI' 

R - ^!^£ P_sPn _ ,P_H^p1. .44. 

l-P-{KB/Ks)+ms/Kf) pI ^ C pI ^^^> 

where Kf,Ks,KB and N are the directly measurable coefficients. This is 
the desired result. 

When ps = or pn = 0, (44) reduces to the set of equations obtained by 
Biot and Willis [8]. Besides, when (3=1 and Kb = N = (the case of bulk 
He II), equation (44) reduces to (18). 



5 Particular solutions 

Equations (28) and (29) describe the propagation of longitudinal and shear 
waves in a superfluid-saturated, homogeneous and isotropic porous media 
in the long wavelength limit when all damping processes are neglected. The 
general solutions of these equations will be discussed elsewhere. Here, we 
consider only the most important particular solutions when 

U = Un- (45) 

Physically this means that the normal fluid component is locked inside a 
porous medium by viscous forces and therefore there is no relative motion 
between the normal fluid component and the porous solid. This case is 
realized in high-porosity aerogels, for example. 
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Because of (45) we added together the first and last equations of (28). 
Then, we obtain: 



V Ph Ph 

= Q^ iipii + 2/>i2 + P22]e + Pu^s) ; 
V2 ([PlQ + R^^]e + R,es) = ^ {pl^e + p^s^s) • (46) 



Ph J dt^ 

Similarly, from (29) and (45) we obtain: 

^^^^ = g^ ([Pll + '^Pl2 + P22]^ + Pu^s) ; 

= Q^iPu^ + P22^s). (47) 

Equations (46) and (47) describe the propagation of longitudinal and 
shear waves respectively, when the normal fluid component is locked inside 
a porous solid. 

It is convenient to introduce the following definitions: 

p' = p + 2^Q + Rn; 
Ph 

Q = Q + Rsn', 

Ph 

p' = p,, + 2p^, + p^,. (48) 

Let us consider equations (46) for the longitudinal waves. The solutions 
of these equations with k wave vector are written in the form 

e = Ciexp[ik{f+Vt)]; 
€s = C2exp[ik{f+Vt)]; (49) 

where V is the velocity of these waves; Ci and C2 are real constants. The 
velocity V is determined by substituting (49) into (46). There are two modes 
of longitudinal waves: the fast and the slow waves, with Vf and Vs velocities 
respectively [1]. The velocities of propagation Vf^s of fast and slow waves 
are given by: 



.2 _ G±[G'-4 {p'p^,, - {pl,f) (P'R, - (Qr)] ^ . 

^^'' Hp'ph-iphf' 



T^2 _ 'J -L [^ - -^ vp P22 - VP12; ; \^ -'t-s - \^ J Ji . /cn\ 
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with 

G = P'pl^ + Rsp'-2Q'p{^; 

The fast wave corresponds to the sohd with locked normal fluid compo- 
nent and the superfluid component moving in phase, and the slow wave 
corresponds to solid with locked normal fluid component and the superfluid 
component moving in opposite phase. Using the formula (50) we examine 
the dependence of Vf^g upon the temperature. In Fig. 6, Fig. 7 and Fig. 8 
fast and slow wave velocities as a function of temperature are presented for 
porosities 0.7, 0.98 and 0.99 respectively (see Section VI for discussion). 
From (46) for fast and slow wave amplitudes we have [1]: 

C2/,s = —Df^sClf^s', (51) 

where 



""'''' Vis p\2-Q" 

Using (15) from (49) and (51) we can write down expression for temperature 
gradient oscillations (recall that e„ = e, in this case): 

AT/^5 = {/STf^s)max exp[i/c(f + Vf^gt) + «vr]; (52) 

where 

Qf~p 

{'^Tf^s)max = ClfyS (1 + Dfys) pr! 

Ph L' 

{ATf^g)max is the temperature oscillations amplitude. The sign of {ATf^g)max 
determines the relative phase of ATj^^ and e (see equations (49) and (52)). 
If {ATf^s)max is positive, ATf^s and e oscillate with opposite phase; and if 
{ATf^g)max is negative, ATj ., and e oscillate in phase. It must be noted that 
{ATs)max is always positive; (ATj)max can be positive or negative. 

Using the formula (52) we examine the dependence of {ATf^g)max upon 
the temperature. In Fig. 6, Fig. 7 and Fig. 8 temperature oscillations 
amplitudes for fast and slow waves as a function of temperature are presented 
for porosities 0.7, 0.98 and 0.99 respectively (see Section VI for discussion). 
We see that {ATf)max is positive in Fig. 6, but is negative in Fig. 7 and 
Fig. 8. 

It is interesting that, in contrast with bulk He II (see Section II), even 
for fast wave e^ 7^ £ni= e) [1] in general. So, fast wave as well as slow wave 
is accompanied by temperature oscillations in general. In the case of bulk 
He II, when /3 = 1 and Kb = N = 0, equation (52) gives that ATf = 0; so 
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temperature oscillations do not accompany first sound, as we have expected 
(see Section II). 

Let us now consider equations (47) for the shear waves. From (47) it is 
easy to see that there is only one type of shear wave [1]. The velocity of 
propagation of these waves Vgh is 



V. 



sh 



N 



p'a-{phy/ip'ph)) 



(53) 



Using the formula (53) we examine the dependence of Vg^ upon the tempera- 
ture. In Fig. 6, Fig. 7 and Fig. 8 the shear wave mode velocity as a function 
of temperature is presented for porosities 0.7, 0.98 and 0.99 respectively (see 
Section VI for discussion). 

Note that, as Pi2 7^ 0, the rotation of the solid with locked normal fluid 
component cj is coupled to the rotation of the superfluid component ujg 
according to the relation 

^s = -4^^; (54) 

P22 

and therefore, ujg ^ [1]. Equation (54) shows that the rotation of the solid 
with locked normal fluid component and the superfluid component are in 
the same direction. 

Equations (50) and (53) give the velocities of propagation of longitudinal 
and shear waves in a superfluid-saturated porous media, when the normal 
fluid component is locked inside the porous solid. If Kf,Ks,KB and N 
will be measured, then one can determine the velocities of propagation of 
fast, slow and shear waves from (50) and (53). It should be noted that 
equations (50) and (53) differ from the equations obtained by Biot [1] and 
contain significant new physics. Indeed, Blot's results do not contain forces 
caused by temperature gradients. In fact, in the case of superfluid-filled 
porous media Biot's results are not valid. Johnson [4] has shown that in the 
rigid frame at low temperatures (T < l.li^) fourth sound is exactly Biot's 
slow wave. However, at low temperatures the hydrodynamical description 
breaks down [3] and comparison to Biot's model should be made at higher 
temperatures [5]. So, we conclude that in all physically interesting cases 
Biot's equations are not applicable to superfluid-filled porous media. 

Examples of physical systems to which equations (50) and (53) can be 
applied are superleak of porous vycor glass and high-porosity aerogels (e.g. 
silica aerogels with porosities /? > 0.9). Superleak of porous vycor glass 
corresponds to a superfluid-solid system in the rigid frame limit (Kf <C 
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Ks,Kb,N). The most interesting physical system, high-porosity aerogel, 
corresponds to a superfluid-solid system in the weak frame limit(i^B, N -C 

6 Discussion and conclusions 

Now, let us discuss the results obtained in the previous sections. In Section 
III and Section IV we derived the general equations (28), (29). These equa- 
tions describe the propagation of longitudinal and shear wave in a superfluid- 
saturated porous media when all damping processes are neglected. Equa- 
tions (28) and (29) depend explicitly on the directly measurable parameters 
of the porous solid and He II via relations (23) and (44). Equations (28), 
(29) and (23), (44) reduce to the equations obtained by Biot [1],[8] when 
Ps = or pn = 0. Besides, when (3 = 1; Kf, = N = these equations de- 
scribe the sound propagation in bulk He II. From what one concludes that 
these equations reduce correctly to Biot's theory [1] and Landau's theory of 
sound propagation in bulk He II [3]. 

In Section V, we considered the case when the normal fluid component 
is locked inside a porous media by viscous forces. It is shown that in this 
case one shear wave and two longitudinal, fast and slow, waves exist. The 
velocities of these waves Vj, Vs, Vgh are obtained. 

Since equations (46) and (47) include pressures caused by thermal gra- 
dients, they are more general in comparison with equations derived by Biot. 
The only case when (46) and (47) reduce to the Biot's equations is p„ = 
(see also Section V). 

Using (50), (52) and (53) we have plotted the examples of the wave 
velocities and temperature oscillations amplitudes for fast and slow waves 
as a function of temperature. In Fig. 6, Fig. 7 and Fig. 8 the plots are 
presented for porosities 0.7, 0.98 and 0.99 respectively. 

The bulk modulus and the mass density of the solid we choose arbitrarily 
as Ks = 4.07 • 10^°iV/m^ and psoi = 4 • IQ^kg/m^. For the bulk modulus of 
the skeletal frame Kb we have used the following relation Kb = {l — fiY'^Ks 
[10]. Poisson's ratio of the skeletal frame is assumed to be 0.2, thus for the 
shear modulus of the solid we have A'^ = (3/4) i^^. When the porosity is not 
small a is given by a = 2 — /? [4] . 

In the present paper the parameters of bulk He II [11] are used for 
numerical estimations. But we note that the thermodynamic parameters of 
Hell in the porous media and the ones of bulk Hell may differ from each 
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other (p„ for example) [3]. 

Figure 6 shows the behaviour of the fast, slow and shear wave velocities as 
a function of temperature for porosity f3 = 0.7. Moreover, the temperature 
oscillations amplitudes for fast and slow waves {ATf)max and {ATs)max 
are presented as well. In this paper we have assumed that the temperature 
oscillations amplitude is about ~ O.IK, so for Ci/ and Cis we choose Cij = 1 
and Cis = 3-10-3. 

The plots presented in Fig. 6 are typical for /3 < 0.9. For these porosities 
Kf <C Ks,Kb,N; so we have a superfluid-solid system in the rigid frame 
limit. In this limit the fast wave and shear wave velocities decrease when 
porosity increases. Besides, the slow wave mode is exactly fourth sound. We 
can show this even analytically. 

In the limit of the rigid frame the expressions (50) for the fast and slow 
wave are greatly simplified. In this case Kf <C Ks, Kb,N and therefore 



v! 



Kb + (4/3)iV 



:i - P)psoi + Pph - {l/a)PPs 



a2 _ 1 (Psj/2 , P«T/2 



Vf = -i'^V.' + '-V,^ ; (55) 

a \Ph Ph } 

here V\ and V2 are given by (21). 

The slow wave corresponds only to the oscillation of the superfluid com- 
ponent (rigid frame) and Vg is the well-known velocity of fourth sound [3] 
renormalized by the tortuous, branching path of pore space [4]; i.e. fourth 
sound is exactly Blot's slow mode. 

In the limit T ^ when there is no the normal fluid component and we 
have {ps/ Ph) = 1) equations (55) are identical to the corresponding expres- 
sions obtained in ref.[4]. 

We see from Fig. 6 that in the fast wave mode {ATf)max > in the 
limit of rigid frame. In the slow wave mode always (ATs)max > (see Fig. 
6, Fig. 7 and Fig. 8). In addition let us note that, in contrast with Vf and 
Vsh, slow wave velocity Vs does not depend on psoi, Kb and N. 

In Fig. 7 and Fig. 8 we present behaviour of the fast, slow and shear 
wave velocities as a function of temperature for porosities /? = 0.98 and 
P = 0.99 respectively. Moreover, the temperature oscillations amplitudes 
for fast and slow waves {ATf)max and {ATs)max are presented as well. For 
Cif and Cis we choose Cij = W~^ and Cis = 3 • 10~^. 

The plots presented in Fig. 7 and Fig. 8 are typical for /? > 0.9. For 
these porosities Kb,N <^ Ks,Kf; so we have a superfluid-solid system in 
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the weak frame limit. In this limit, in contrast with (3 < 0.9 (rigid frame), 
the slow wave velocity Vs depends on psoi, Kb and N. 

Fig. 7 shows that Vs is less compared with fourth sound velocity, but 
still is larger than second sound velocity. 

In Fig. 8 we see that the behaviour of Vg near Tx is almost the same as in 
Fig. 7. But at low temperatures Vg becomes considerably less than second 
sound velocity. For porosities f3 > 0.99 slow wave velocity Vs increases 
smoothly when f3 increases; and when porosity reaches (3 = 1 the slow wave 
mode is exactly second sound. 

Fast wave velocity Vf decreases when the porosity increases for /3 < 0.95 
(it becomes even less than first sound velocity); and Vj increases when the 
porosity increases for (3 > 0.95. In Fig. 7 and Fig. 8 Vf is less than first 
sound velocity. For (3 = 1 fast wave is exactly first sound. The shear wave 
mode velocity Vgh decreases when the porosity increases. Note that in Fig. 
7 and Fig. 8 Vsh is less than Vs. For (3=1 the shear wave velocity is zero. 

From Fig. 7 and Fig. 8 we can see that in the fast wave mode {ATf)max < 
in the limit of weak frame. When porosity is exactly (3 = 1 (the case of 
bulk He II), the temperature oscillations amplitude in the fast wave mode 
is zero (also see Section V). 

The limit of weak frame is the most interesting and important application 
of our theory. Indeed, in recent years the investigation of the behaviour of 
superfluids in aerogels has attached increasing attention (for example, see 
[12], [13], [14]). Aerogels are highly open porous solids with porosities (3 > 0.9 
(e.g. silica aerogels). The mean pore diameter in aerogels is about 20nm, 
so the normal fluid component is clamped in the frame by viscous forces. 
Besides, high-porosity aerogels have very small bulk and shear modulus and 
are highly compliant. That is why in the superfluid-filled aerogels pressures 
caused by thermal gradients play a crucial role, in contrast to other porous 
solids. Thus, one should describe wave processes in high-porosity aerogels 
using equations (46) and (47). 

We point out that the theory of sound propagation in aerogels has been 
developed in ref.[12]. These authors modified the conventional two fluid hy- 
drodynamic equations and derived an expression (see equation (8) in ref.[12]) 
which gives the longitudinal wave velocity in superfluid-filled aerogels, but 
does not give the shear wave velocity. It must be noted that equation (8) in 
ref.[12] is incorrect. To show this we remark the following. 

The longitudinal wave velocity in any solid (in aerogel in the case of 
ref.[12]) is given by [{Ka + i4/3)Na) / Pa]^^^; [15]. Here Ka, Na and pa are 
the bulk modulus, the shear modulus and the mass density of the solid. 
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respectively. In [12] the term {A/2>)Na/ pa is omitted (see the equation (5) 
in ref.[12]). This is justified if Na/Ka «C 1. In the case of aerogels Poisson's 
ratio equals ^ 0.2 and is independent of density [16], [10]. This means that 
Na/Ka ~ 1 and the term {4:/3)Na/pa in the longitudinal wave velocity is 
not negligible. 

Besides, in the limit pa ^ oo equation (8) in ref.[12] suggests that in the 
superfluid-filled porous media one of the modes of sound propagates with 
velocity Ca, where Ca is given by equation (5) in ref.[12], which is incorrect 
(see [15]). 

In conclusion, we have derived the general equations which describe the 
propagation of longitudinal and shear waves in a superfluid-filled, homo- 
geneous and isotropic porous media in the long wavelength limit when all 
damping processes are neglected. The most important case when the normal 
fluid component is locked inside a porous media by viscous forces is inves- 
tigated in detail. It is shown that, in this case, one type of shear wave and 
two types of longitudinal, fast and slow wave, propagation are possible. The 
velocities of these waves are obtained. Our results reduce to the well-known 
equations in the appropriate limits. We have also demonstrated that the 
formalism developed by Biot may be applied to the wave processes in He 
II. It is shown that fourth sound, as well as second sound in bulk He II is 
exactly Blot's slow wave. In addition, we remark that the results derived 
in this paper may be applied to the porous media saturated with superfluid 
^He as well. 
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FIGURE CAPTIONS : 



Fig. 1. The "unjacketed compressibility test" for bulk He II. A pressure 
p is applied to a volume of He II in a vessel. In this case none of the fluid 
components is allowed to escape from vessel. 



Fig. 2. The "second jacketed compressibility test" for bulk He II. A 
pressure p' is applied to He II in a vessel. The applied external pressure 
squeeze out some of the superfluid component through the superleak into 
the reservoir of He II. As a consequence temperature difference between the 
fluid in the vessel and the reservoir of He II exists. 



Fig. 3. The "unjacketed" compressibility test. A pressure p is applied to 
a sample of the superfluid-filled porous solid. In this case none of the fluid 
components is allowed to escape into the reservoir of He II. 



Fig. 4. The first "jacketed" compressibility test. A pressure p is applied 
to a sample of the superfluid-filled porous solid. The applied external pres- 
sure squeeze out some of the fluid through the tube into the reservoir of He 
II. In this case both of the fluid components are allowed to escape into the 
reservoir of He II. Thus, no temperature difference between the superfluid 
in the porous solid and the reservoir of He II exists. 



Fig. 5. The second "jacketed" compressibility test. A pressure p is ap- 
plied to a sample of the superfluid-filled porous solid. The applied external 
pressure squeeze out some of the superfluid component through the super- 
leak into the reservoir of He II. In this case only the superfluid component is 
allowed to escape into the reservoir of He II. As a consequence temperature 
difference between the superfluid in the porous solid and the reservoir of He 
II exists. 



25 



Fig. 6. The velocities and the temperature oscillations amplitudes as a 
function of temperature for fast and slow waves; and the shear wave velocity 
as a function of temperature. (3 = 0.7, Cif = 1 and Cis = 0.003; Ks = 4.07- 
W^^N/m^ and psoi = 4 • lO^kg/m^; Kb = {I - ^f'^Kg and A^ = (3/4)Kb; 
a = 2 — (3. Here the parameters of bulk He H are used for numerical 
estimations. 



Fig. 7. The velocities and the temperature oscillations amplitudes as a 
function of temperature for fast and slow waves; and the shear wave velocity 
as a function of temperature. (3 = 0.98, Cif = 0.1 and Cis = 0.3; Ks = 4.07- 
W^^N/m^ and p^oi = 4 • lO^kg/m^; Kb = {I - Pf'^Kg and A^ = (3/4)Kij; 
a = 2 — (3. Here the parameters of bulk He H are used for numerical 
estimations. 



Fig. 8. The velocities and the temperature oscillations amplitudes as a 
function of temperature for fast and slow waves; and the shear wave velocity 
as a function of temperature. (3 = 0.99, Cif = 0.1 and Cis = 0.3; Ks = 4.07- 
lO^OiV/m^ and Psoi = 4 • lO^^/m^; Kb = {1- l3f-^Ks and N = (3/4)Kb; 
a = 2 — (3. Here the parameters of bulk He II are used for numerical 
estimations. 
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